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The behavior of perturbation in scalar field dark energy and its consequent effect on the cold dark 
matter (CDM) power spectrum is well understood to be governed by the equation of state (EOS) 
parameter and the effective speed of sound (ESS) of dark energy. In this paper, we investigate 
whether dark energy models whose ESS are epoch dependent leaves any distinct imprints on the 
large scale CDM power spectrum. In particular, we compare the cases where the ESS is decreasing 
with time with those where it increases. The CDM power spectrum is found to be generically 
suppressed in these cases as compared to the ACDM model. The degree of suppression at different 
length scales can, in principle, refiect the evolving nature of the ESS of dark energy. However, we 
find that the effect on the CDM power spectrum in cases where the ESS of dark energy is evolving 
with constant EOS parameter is significantly smaller as compared to the situation where ESS is 
constant whereas EOS parameter is evolving. Further, it is also shown that the effect of different 
, evolution of ESS for a given evolution of EOS parameter of dark energy on the CDM power spectrum 

'^'n , is significant only at the intermediate scales (around k ~ O.Olh/Mpc). At scales much smaller and 

' larger than the Hubble radius, it is the evolution of EOS parameter of dark energy which governs 

I the degree of suppression of CDM power spectrum with respect to the ACDM model. 

PACS numbers: 98.80.-k, 95.35+x, 98.65.Dx 
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p I. INTRODUCTION 

O . 

$-H , It is now evident from numerous cosmological observations that our Universe is undergoing accelerated expansion 
C/3 ■ at the present epoch [1-5]. This profound discovery triggered a flood of investigations both at the theoretical and 
I experimental front to understand the reason for this accelerated expansion. The theoretical explanation for this 
observed late time accelerated expansion can broadly be classified into two categories. The first explanation is based 
on the assumption that gravity as described by the Einstein's general theory of relativity (GTR) is valid on the 
J> . cosmological scale. This assumption then forces us to invoke some exotic matter known as dark energy to drive the 
On ■ late time accelerated expansion of the universe. If not the dark energy, then the only other way of explaining this 
' accelerated expansion is to modify the theory of gravity from the standard GTR. In the present work we shall only 
consider the first case, viz., the one where the accelerated expansion is driven by dark energy. 

Numerous models of dark energy have been investigated in the literature, these include quintessence, tachyon, 
. k-essence, Chaplygin, gas etc. (for details, see Refs.[6-9] and the references therein). It is now a big challenge to 
' observationally test the predictions of these models. One of the ways of distinguishing various models of dark energy 
, is to investigate how cold dark matter (CDM) clusters on cosmological scales. In this paper we shall investigate how 
. . ' a general class of scalar field models of dark energy influences the large scale CDM power spectrum. 
K*" It is well known that scalar field models of dark energy can lead to degenerate evolution of the scale factor [10, 11]. 

Therefore, the evolution of scale factor is insufficient to distinguish scalar field models of dark energy. Hence, it is 
'^^j ! also necessary to take into account the consequence of perturbation in dark energy on cosmological observable such 
d . as CDM matter power spectrum. Although, the evolution of scale factor and the growth of CDM perturbations does 
not uniquely characterize the nature of scalar field lagrangian of the dark energy (see Ref.[12] for details), however, 
such studies can be useful in at least ruling out some of the models. 

Evolution of perturbations in dark energy and how it infiuences CDM power spectrum and the integrated Sachs- 
Wolfe (ISW) effect in the cosmic microwave background (CMB) have been investigated in the literature [13-36]. It 
is well understood from these studies that CDM power spectrum is suppressed in a generic scalar field model of dark 
energy as compared to the same in ACDM model. The degree of suppression depends on how scale factor as well as 
perturbations in dark energy evolves. It is important to take into account the role played by the perturbations in 
dark energy for the system of equations to be consistent with Eistein's equation [23, 26]. 
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In general, any minimally coupled scalar field model of dark energy can be characterized by the behavior of the 
following two parameter: (1) equation of state (EOS) w, and (2) the effective speed of sound (ESS) c,.- Most studies 
on the influence of perturbation in dark energy involves the case where Cg and w are constants (see for instance, 
Refs.[30, 35]). Perturbations in dark energy models where ESS parameter Cg is epoch dependent have also been 
investigated [16, 33]. However, the evolution of Cg was fixed by the specific model of dark energy under investigation. 
Our aim in this paper is to parameterize the evolution of Ce as a function of scale factor and to study how perturbations 
in such dark energy models influences the CDM power spectrum. 

In this paper, we shall consider a class of k-essence dark energy models with lagrangian of the form £(X, (j)) = 
F{X) — V{(f)), where X = (1/2) 9^ is the kinetic term. We shall first present a method of reconstructing the form 
of F{X) and V{(l)) from a given evolution of EOS parameter w(a) and ESS parameter Ce(a) with scale factor. Using 
this method, we shall reconstruct few models of dark energy starting from a desired evolution of w{a) and Ce(a). The 
evolution of perturbation in dark energy is investigated for the following cases (i) constant w and Ce, (ii) either of w 
and Ce is evolving, and (iii) both w and Cg evolving. The growth of large scale CDM perturbation is compared in 
these models. It is shown that in models where w is constant whereas Cg is evolving, the effect of different functional 
form of Ce(a) on the CDM power spectrum is significantly smaller as compared to the case where Ce is constant and 
w is evolving. 

This paper is organized as follows. In the next section, we discuss the necessary equations required for studying 
the background evolution in a system of CDM and k-essence dark energy. Section III deals with the evolution of 
perturbations in the longitudinal gaiige in a system consisting of CDM matter and scalar field dark energy. In section 
IV, we develop a method of reconstructing the Lagrangian density of the scalar field dark energy from given the 
evolution of equation of state parameter and sound speed of dark energy perturbations. The influence of dark energy 
perturbations on the CDM power spectrum in such reconstructed models is studied in section V. Finally in Sec. VI, 
we discuss in detail the effect of ESS parameter at different length scales in CDM power spectrum. Summary and 
conclusions are given in section VI. 



II. BACKGROUND EVOLUTION 



We consider a spatially flat homogenous and isotropic Friedmann-Robertson- Walker (FRW) universe with line 
element given by 

rf.s^ = dt^ - (t) [dx^ + dy^ + dz"^] , (1) 

where a{t) is the scale factor. The late time evolution of the universe is primarily governed by the properties of dark 
matter and dark energy. Hence, the late time evolution of the scale factor is determined by the following Friedmann 
equation 

[Pm{a) + Pde{a)\ , (2) 

where pm and pd are the energy densities of dark matter and dark energy respectively. If the dark matter and dark 
energy are not coupled to each other, then each component individually satisfies the following conservation equation 

p = -3H{p + p). (3) 

Since the dark matter is non-relativistic (p ~ 0), the conservation equation [Eq.(3)] implies that 

Pm.=Pmoa~^, (4) 

where pmo is the density of the dark matter at the present epoch. From the conservation equation (Eq.(3)), it follows 
that the dark energy density evolves as 

Pde{a) = pdeo exp ( -3 / n + ^(") j ^ (5) 



where pdeO is the dark energy density at the present epoch and u;(a) is the equation of state parameter of the dark 
energy defined as 

w{a) = ^. (6) 

Pde 
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For a given w{a), Eqs.(2), (4), and (5) form a close set of equations determining the evolution of scale factor. 
However, the evolution of the equation of state parameter w(a) with scale factor can only be determined by the 
underlying nature of dark energy. 

We assume that the dark energy is a minimally coupled scalar field with the action given by 



S[cj>] 



I 



-g jC{X,cj>)d''x, 



where C{X, </>) is the lagrangian density of the scalar field which is a function of scalar field 
X, which is defined as 

X = {l/2)d^<f>d''<p. 

The stress energy-momentum tensor associated with the scalar field is given by 



(7) 

and the kinetic term 



dC 
dX 



{d^(l>d''(t>)- g'"' C. 



In the background space-time with the FRW line element, the scalar field (j) can only be a function of time < 
Consequently the stress energy-momentum tensor will be diagonal T^" = diag(p^, — p<^, — p<^, —P(p), with 

Therefore, the equation of state parameter of dark energy defined in Eq.(6) can be expressed as 



w ■ 



2{^)X-C{X,<l>y 



(8) 

' = m- 

(9) 
(10) 



(11) 



The above equation [Eq.(ll)] implies that the evolution of EOS parameter w{a) can be determined once the evolution 
of the scalar field (t>{t) on the FRW background is known. The scalar field action (7) implies that the background field 
<j){t) satisfies the following equation of motion 



dx) 



2X 



\dX-^ 



3H 



\dX 



dXdcj) 



= 0. 



(12) 



Given a Lagrangian density £(X, 0), Eqs.(2) and (12) form a closed set of equations for determining the evolution 
a(t) and (/'(t). Consequently, the evolution of w(a) can be determined using Eq.(ll). However, it is also possible to 
reconstruct the form of Lagrangian density starting from from a given evolution of the EOS parameter w(a). We 
shall adopt this method of reconstruction in section HI for designing few models of dark energy which leads to some 
desired evolution of w{a). 



III. EVOLUTION OF SCALAR PERTURBATION 



In order to study the evolution of perturbations in CDM and dark energy we consider the following perturbed FRW 
metric in the longitudinal gauge [37-39], 

ds^ = (1 + 2$) dt^ - a^it) (1 - 2$) [dx^ + dy^ + dz^] , (13) 

where $ is the variable describing scalar metric perturbations and it is also known as the Bardcen potential [37] . The 
energy momentum tensor for the matter content of the universe such as CDM and scalar field dark energy can be 
expressed as 

T'', = {p + p)u^'u,-p6^',, (14) 

where p is the energy density, p is the pressure and is the four velocity field. The perturbations in the energy 
density, pressure and the four velocity field are defined in the following way 

pit,x) = pit) + Spit, x), (15) 
p{t,x) = Pit) + Spit, x), (16) 
u^" = u^' + Su^', (17) 
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where p{t) and p{t) arc the energy density and pressure in the background FRW hue clement, respectively, and 



= [1,0,0,0] is the background four velocity. Since 



1, it turns out that 6uo 



spatial part of the four velocity field Su'' can be expressed as a gradient of a scalar 



-6u' = Further, the 



(18) 



From Eqs.(14) to (18), it follows that the perturbations in the energy momentum tensor can be written in terms of 
6p, dp and u as 

ST°o = Sp, 

5T\ = {p + p)ui 

The three variables Sp, dp and u describes the scalar degree of perturbations in the matter sector for both perfect 
fluid and scalar field. In the case of pressureless matter {p = 5p„i = 0), the perturbations are described by 5pm and 
Um- For the scalar field dark energy, the perturbation in the scalar field <p is defined as 



(l){x,t) = (l){t) + 6(j){x,t), 



(19) 



where (j){t) is the value of scalar field on the background FRW spacetime. Substituting Eq.(19) in the stress-energy 
momentum tensor of scalar field tensor defined in Eq.(8) and subtracting the background and p^ we get 



f dC 



'd4> 



dL 




(20) 
(21) 

(22) 



In the system consisting of CDM and dark energy, the linearized Einstein's equation SC^^^ = SwGST^^, which 
relates the matter perturbations to the metric perturbations, leads to the following equations 



-AttG {pm{a)5m + Pd{a)Sd) ■. 



^ + H^ = -A-kGo^ {Pm{a)Um + {Pd + Pd)Ud) , 
$ + 4i?6 + (2H + SF^) $ = 47rG Spd, 

where Sm and 84 are the fractional density perturbation^ in CDM matter and dark energy, respectively, defined as 

r _ Spm 

Om = — 



(23) 

(24) 
(25) 



5d ^ ^. 

Pd 

It should be understood from the above Einstein's equations [Eqs.(23)-(25)] that each of the variables describing the 

perturbations such as 5pm, Spd, etc. are in fact the amplitude in the fourier space for a given fourier mode k. 
The covariant conservation equation T^^ ^.^ = implies that for the pressureless matter such as CDM 



5r, 



= -2HUrr,- 



^ Hereafter, we will be using the subscript 'd' instead of ij> in the variables describing perturbations in dark energy 
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where overdot denotes derivative with respect to cosmic time. In the case of scalar field dark energy the perturbations 
in pressure 6pd is related to 6pd and as [12, 40, 41] 

5pd = clSpd - 3H {pd + Pd) a^Ud [cl - cl] , (26) 

where is the square of the adiabatic speed of sound defined as 

2 _ Pd W 

c„ = — = w — -rrrrr-. (27) 

and Cg is square of the effective speed of sound of scalar field given by [42, 43] 

(M.\ 

c2 - \ml (28) 

\dx) \dX^) 

The ESS parameter Cg is in fact the ratio of dark energy pressure perturbation to the energy density perturbation in 
the comoving gauge or rest frame of dark energy [18, 44]. In the case of scalar field dark energy this gauge coincides 
with the uniform field gauge. 

For the dark energy density perturbation, the covariant conservation equation {T'^v;ij, = 0) together with Eq.(26) 
leads to the following equations 

5d = {1 + w) k^Ud + 3H{w- cl) Sd + 9H^ (1 + w) [cl - cl] a'^Ud + 3 (1 + «;) (29) 

iid = -H{2- 3cl) Ud - ,ff' . - * . (30) 

We now introduce the following dimensionless variable cTd defined as 

ad = 3H{l + w)a^Ud (31) 

For studying the evolution of perturbations in a system of CDM matter and dark energy, it is evident that Eqs.(25), 
(29) and (30) forms a closed set of equations. In terms of the variables Sd and ad, these equations can be re-expressed 
as 

S'd = (^)(---e) + ^(c^-c^) + ^+3(l + H^' (34) 

where prime denotes the derivative with respect to the scale factor a. Once the solutions $(o, k) and 5d{a, k) are 
determined from the above three equations, the evolution of perturbation in CDM can be evaluated from the time-time 
component of the Einstein's equation and this is given by 

« = - (a^) + (siS?) *) - iM) '''' 

where Q.m and Q.d are the dimensionless density parameter, defined as 



IV. RECONSTRUCTING DARK ENERGY MODELS 



It is evident from Eqs.(32)-(34) that the evolution of perturbations in the matter and dark energy are governed by 
the behavior of dark energy speed of sound Cg(a) and equation of state parameter w{a). As mentioned earlier the 
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evolution of Cg(a) and w(a), in general, depends on the underlying nature of the lagrangian density of the scalar field 
describing the dark energy. However, one could also reconstruct the form of the lagrangian density of the scalar field 
dark energy from a given evolution of ESS and EOS parameter. In this section we will present such a method of 
reconstruction. The evolution of perturbation in such reconstructed models will be discussed in the following section. 

In this paper we restrict our analysis to the following class of scalar field dark energy models with the lagrangian 
density given by [45]: 



C{X,<I>) = F{X)-V{<I>). 



(36) 



This is a natural generalization of the quintessence model where F{X) = X. The form of the kinetic function F{X) 
and potential V{(l)) can be reconstructed from a given evolution of w{a) and Cg(a). It is important to note that such 
class of dark energy models allows evolution in Cg(a) even when w is constant. This is unlike the tachyon model 
{C{X,(j)) = —V{(p)\/1 — 2X) where Cg(a) = —w{a) and consequently, a constant equation of state always necessarily 
implies that cf, is also constant. 

For the model described by the lagrangian density (36), the field equation for the (j) given by Eq.(12) becomes 



\Fx 



Fx \ 

- 2XFxx J 



Fx + 2XFxx 



= 0. 



The ESS parameter Cg for this model reads 



^x 



Fx + 2XFxx 



Further, the EOS parameter w{a) can be expressed as 



w = 



2XFxx - jC 

Using Eqs.(38) and (39), the equation of motion for the scalar field (p can be re-expressed as 



X 



6H + 



w 



V 



Pd{a) 



(37) 



(38) 



(39) 



(40) 



where X = ^cf)^ \s the kinetic term for the backgrormd field 



From the definition of adiabatic speed of sound (see Eq.(27), it follows that 



V -FxX 
GHXFx ' 



(41) 



which on rearranging implies that 



X 
X 



V 



1 + w. 



Pdia) 



(42) 



The two equations [Eq.(40) and (42)] are the evolution equations for the scalar field, with the first one being the 
field equation derived from the scalar field action and second one follows from the definition of the adiabatic speed of 
sound. Equating the two equations we get 



dV 
da 



1 + d 



w[a] 



exp 



1 + w{a) 



da 



(43) 



X{a) = Xoexp 



■I- 



da 
a 



(44) 



where Xq is the constant of integration. Since, X = hfj)^, the evolution of scalar field with the scale factor is given by 



(45) 
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FIG. 1: The percentage suppression in matter power spectrum in dark energy models in comparison to ACDM model is plotted 
for a constant equation of state w and for different values of c^. The left panel is for w = —0.9 and right panel shows suppression 
for w = —0.8. It can be seen that suppression is more for dark energy models with w = —0.8 with respect to w = —0.9. The 
power spectrums shown in this figure corresponds to its value at the present epoch. 



Given any functional form of Cg(a) and w(a), Eq.(44) can be integrated to obtain X(a). Consequently, the evolution 
of scalar field with the scale factor (p{a) can be determined from Eq.(45). Further, integration of Eq.(43) would give 
potential V{a) as a function of scale factor. By eliminating the scale factor 'a' from the two functions V{a) and (f>{a), 
the form of the potential V{(l)) can be determined. 

The form of the kinetic function F{X) can be reconstructed in the following way. From the definition of the EOS 
parameter (see Eq.(39)), it follows that 



dF 
d\nX 



1 + w{a) 



Pd{a). 



(46) 



Using Eq.(44), the above equation can be re-expressed as 



dF 
da 



1 



3(l + ii;) 



Prfo exp 



1 + w{a) 



da 



(47) 



From the solution of the above equation, it is possible to determine how the kinetic function F{a) evolves with scale 
factor, for a given evolution of w{a) and Cg. Further, from the two solutions: X{a) obtained in Eq.(44) and F{a) 
from Eq.(47), one can eliminate the scale factor to reconstruct the form of the kinetic function F{X). It may not 
always be possible to reconstruct the exact analytical form of F{X) and V{(l)), however, one can always numerically 
determine the functional form of these functions. As our primary interest in this paper is to investigate how dark 
energy perturbations influences the CDM power spectrum, we will not be discussing the numerical form of F[X) and 
V{(j)) in each and every case. The main purpose of this section was to emphasize the fact that the Lagrangian density 
of the form £(X, 0) = F{X) — V{4)) indeed allows a class of solutions with evolution in Cg(a) and w{a). 



V. SUPPRESSION OF CDM POWER SPECTRUM 



Our aim in this paper is to investigate how perturbations in the dark energy influence the CDM matter power 
spectrum. In particular, we are interested in models of dark energy whose speed of sound is epoch dependent. 
In the preceding section, it was illustrated that it is possible to reconstruct the Lagrangian density of the form 
C{X, (f) = F{X) — V{(j)) from a given evolution of the EOS and ESS parameters. Hence, in principle, it is always 
possible to associate F{X) and V{(f>) for a given evolution of Cg(a) and w{a). This allowed us to choose a parameterized 
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FIG. 2: Plot shows the ratio of the percentage suppression in matter power spectrum A defined in Eq. (48) for w = —0.9 to 
w = —0.8, for different values of effective speed of sound. 

form of Cg(a) and w{a) for which the evolution of perturbations can be investigated. But, before proceeding to the 
case where Cg(a) is epoch dependent, we will first discuss the case where it is constant. 

Case I: Constant c? and w 



Considering the class of scalar field dark energy models with Lagrangian density of the form C{X, i 
if both w and Cg(a) are constants, it follows from Eqs.(44) and (47) that X{a) oc and F{a) (x a^' 
n2 are constants. Therefore, when both w and c; 
would become [45] 



:F(X)-F(0), 

where ni and 

:(a) are constants, the functional form of the kinetic function F{X) 



F{X) ^X" 



where a is a constant given by 



2c? 



It should be noted that for models with F{X) = X", the ESS parameter Cg is always constant irrespective of the 
functional form of the potential V{(j)) and w{a). It is the form of the potential which would lead to a constant value 
for the equation of state parameter. The form of the potential can be numerically determined from Eqs.(43) and (45). 

In order to determine how CDM perturbation grows when both Cg and w are constants, we numerically solve the 
closed set of equations (32) to (34), which takes into account the role played by the perturbations in dark energy. 
Once the evolution of $(a), 5d{a) and (Td(a) is known, the growth of CDM perturbations can be determined using 
equation (35). For quantifying the difference in the growth of CDM perturbation in dark energy models with respect 
to ACDM model, we numerically evaluate the following quantity 



A% 



QCDM 



(fc) 



,(fc) 



P. 



X 100 



(48) 



where Pqchm (^) is the CDM power spectrum at the present epoch in dark energy models whereas Pj^com (^) the 
corresponding power spectrum in ACDM model. This quantity for the case of constant and w is plotted in Fig. 1. 
It is clear from this figure that A is negative at all length scales of perturbations. Therefore, CDM power spectrum 
is suppressed in dark energy models as compared to the same in ACDM model. 
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FIG. 3: Figure showing suppression in matter power spectrum compared to ACDM model for the case when w is evolving, but 
Ce is constant. In the left panel, w is initially —1 and asymptotically approaches —0.8 whereas in the right panel w is initially 
—0.8 and asymptotically approaches —1. It can also be verified that in both the panel w at the present epoch is —0.9. 



The degree of suppression in CDM power spectrum depends on the value of w and Cg and on the length scale of 
perturbation. In Fig. 1, we have plotted the percentage suppression A% for w =^ —0.9 and w — —0.8 and for different 
value of constant Cg. It follows from this figure that there is more suppression if w — —0.8 than if = —0.9, indicating 
the fact that the degree of suppression is in a way proportional to the deviation of w from —1. When comparing 
percentage suppression for w — —0.8 and w = —0.9 in Fig. 1, it may appear that they differs just by a factor which is 
independent of Cg. However, this is not true as illustrated by Fig. 2, where we have plotted the ratio of A (defined in 
Eq.(48)) for w — —0.8 and w — —0.9. This ratio does depend on the value of clarifying the fact that the percentage 
suppression A is not merely a product of two functions with one depending on w and the other on Cg . 

Further, for any value of w and Cg, the percentage suppression decreases with increasing length scale of perturbations. 
This statement is true for the case of constant w, however, in general, it depends on the functional form of the equation 
of state parameter w{a). 

It can also be observed from Fig. 1 that the effect of different values of Cg (assumed to be constant) for a given value 
of w on the CDM power spectrum is insignificant on scales much smaller than the Hubble radius. This statement 
is also true for scales much larger than the Hubble radius. In fact, it can be be seen from Fig. 1 that only at the 
intermediate scales, say at fc = O.Olh/Mpc, the effect of different Cg is more pronounced. At these scales, it turns out 
that the percentage suppression increases with the increase in the value of Cg . The reason for the fact that the effect 
of different values of Cg is significant only at the intermediate scales will be discussed in detail in the next section 
(Sec. VI). 



Case II: Constant and evolving w 

When Cg is constant, it is in general possible that the equation of state parameter is epoch dependent. We will now 
consider such a case. As mentioned earlier, for models with F{X) — X", Cg is always a constant. In fact, it is the 
form of the potential V{(j)) in these models which determines the evolution of the EOS parameter. We consider the 
following parametrization for the evolution of w{a), 

w(a) = Wo + wi (^Y"^) (49) 

where Wq, Wi and n are constants. The constant wq correspond to the initial value of the equation of state parameter 
(at a = 0). Asymptotically the equation of state parameter approaches the value wq + wi. The parameter 'n' describes 
the rate at which w evolves from its initial value wq to its asymptotic value wq + wi. The value of w{a) at the present 



10 




FIG. 4: Figure showing suppression in matter power spectrum for w{a) given by Eq.(49) with n = 10. In the left panel the 
equation of state parameter is initially —1 and stays close to it most of the time whereas in the right panel it evolves from 
—0.8. If one compares the suppression at two length scales say at = O.lh/Mpc and at = O.OOOlh/Mpc, it turns out that 
suppression is more at larger scale compared to small scales in the left panel of this figure. The situation is exactly reverse in 
the plot on the right panel. This illustrates the fact that the behavior of the suppression at these scales precisely depends on 
the functional form of w{a). 



epoch would be wq + {wi/2"'). Hence, greater the value of n slovifer will the evolution w from wq to its asymptotic 
value wq + wi. 

The parametrization Eq.(49) is different from the one generally investigated in the literature, viz. the Chevallier- 
Polarski-Linder (CPL) parametrization [46, 47] where w(z) ~ wq + w^ (^TTz^ ■ However, at low red shifts the functional 
form of both the parametrization converges. Although, the CPL parametrization is very well valid for the range 
of red shifts which is observationally relevant, it does not restrict the value of w to less than unity. In fact in 
CPL parametrization, the equation of state parameter diverges asymptotically. This problem does not arise in the 
parametrization of w{a) introduced in Eq.(49). 

For constant value of Cg, we assume that w{a) evolves with scale factor with the functional form given in Eq.(49). 
In order to study the effect on matter power spectrum by the dark energy perturbation, we consider two specific 
behaviors of w{a): one in which w{a) is initially -1 and approaches a constant value {w > —1) asymptotically. The 
other case of interest is where w > —1 initially, but asymptotically approaches -1. For the case where w is initially -1 
we choose wq = —1 and wi to be 0.2, such that asymptotic value oi w = —0.8. The range of w between -1 and -0.8 is 
well within the observable constraint on the value of w. The second case of interest where w decreases with time and 
finally approaches -1 would require Wq + Wi — —1. We assume that initially w — —0.8, therefore Wi would be -0.2. 

The suppression of CDM power spectrum as compared to the ACDM in these two cases is shown in Fig. 3 for value 
n = 1 in Eq.(49). As in the case of constant of w and Cg, here also we see that there is a suppression at all scales. 
The suppression is more in the case when w = —0.8 initially (around 8 percent) than the case where w is initially -1. 
Again it can be seen from Fig. 3 the effect of different values Cg on the suppression is more pronounced at intermediate 
values of comoving wavenumber k. 

In the above cases we have chosen the parameter n in Eq.(49) to be unity. If we increase n to 10, then the rate at 
which EOS parameter switches its value from initial value to the asymptotic value wq + wi slows down dramatically. 
This case is shown in the Fig. 4. When n — 10, the equation of state parameter for most the time stays near its initial 
value. Therefore, if wo = ~1, then it nearly mimic ACDM model for most of evolution, but deviates from it at the 
present epoch. The peculiar behavior of the suppression in the CDM power spectrum in this case is shown in the left 
panel in Fig. 4. Also we note from this plot that the suppression is very minimal in this case and contrary to other 
cases the suppression is maximum at small fc's than large values of k. This peculiar behavior is due to the choice of 
the function w{a). 

The right panel of Fig. 4 shows the other case where w is initially -0.8 and approaches -1, the percentage suppression 
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FIG. 5: The behavior of suppression in matter power spectrum as compared to ACDM model for the case when EOS parameter 
w is constant but ESS parameter is evolving with scale factor is shown in this figure. The left panel shows the case when 
w — —0.8 whereas the right panel show for w = —0.9. The suppression in matter power spectrum for the case when 
asymptotically decays to zero from its initial value of unity is compared with the reverse case where it increases to unity from 
zero. For comparison = 0.5 is also plotted. Note that in all these cases the value of Ce at the present epoch is 0.5 



in this case is high for the obvious reason that w for most of the time stays close to —0.8. Also it should be noted 
that suppression in this case is larger at smaller length scales unlike the case in the left panel. This illustrates the 
fact that the percentage suppression at both smaller and larger scales very much depends on the functional form of 
the equation of state parameter w{a). It should also be noted that the effect of different values of Cg can be seen more 
profoundly only at the intermediate scales. 



Case III: constant w and evolving Cg 

Let us now consider the case of dark energy models for which EOS parameter w is constant but evolves with 
time. As mentioned earlier, the class of dark energy models described by the lagrangian density (36) does allow to 
evolve with time when w is constant. It should be noted that this is not in general valid for all classes of scalar field 
DE models. For example in tachyon model of DE with Ci — V^(0)\/1 — 2X, a constant w always implies that Cg is 
also constant. 

We consider the following parametrization for the evolution of Cg with the scale factor 

=-^<* + '-(tIi)" <'^'" 

This is similar to the parametrization of w discussed in the preceding section (see Eq.(49)). The interpretation of the 
constants cq and ci in Eq.(50) is as follows. The constant cq corresponds to the initial value of the ESS parameter 
Cg, whereas cq + ci corresponds to its asymptotic value. The parameter n determines the rate at which the ESS 
parameter evolves from cq to its asymptotic value cq + ci . 

We are interested in comparing the suppression in CDM power spectrum with respect to ACDM in cases where Cg 
decreases with time with those where it increases. For the increasing case we choose cq = and ci = 1 so that Cg 
evolves from to 1 asymptotically. Similarly for the decreases case we choose cq = 1 and ci = — 1 so that Cg evolves 
from 1 to 0. The suppression of CDM power spectrum for these cases with n = 1 in Eq.(50) is shown in Fig. 5. This 
choice of these set of parameters ensures that the value Cg at the present epoch for these two cases is the same and 
is given by 0.5. Hence, for comparison the case of constant Cg with Cg — 0.5 is also plotted in the same figure. It 
is interesting to note that the three curves in Fig. 5 (both left and the right panel) almost overlap with each other 
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FIG. 6: Comparison of suppression of matter power spectrum for two different evolution of Cg(a) with scale factor is sfiown in 
tliis Fig. Here we have taken co — 10"^ and n = 10. For both these curves it can be verified that the present value of is 
approximately 10~^. 

indicating the fact that the effect of different evolution of Cg (a) for a constant w is insignificant on the suppression of 
CDM power spectrum. This may not be true for any evolution of Cg with the scale factor. 

In order to check the generality of this result we consider two more different evolution of Cg. The first one is the 
same function introduced in Eq.(50) with the following values of the parameter cq = 10^^, ci = 1 and n = 10. In 
this case Cg(a) evolves from cq = 10~^ initially to 10~^ at the present epoch. The second functional form of Cg(a) we 
consider is Cg(a) = 1/(1 + a)" with n — 10. Here, in this case Cg(a) evolves from its initial value of unity to 10~^ at 
the present epoch. The suppression of CDM power spectrum in these two cases is shown in Fig. 6 It is evident from 
this figure that the two curves are not as close to each other as compared to the three curves in Fig. 5. The reason for 
this is that in these two curves order of magnitude change in Cg from z = 1000 to present epoch is much larger than 
the cases presented in Fig. 5. However, it should be noted that, as in the previous cases of constant Cg the effect of 
different choice of Cg(a) on CDM power spectrum is not evident at large and low values of k. 

Case IV: Evolving w and 

Until now we have considered the cases where both w and Cg are constants and the case where either one of them 
is epoch dependent. For the purpose of completeness, we now consider the case where both w and Cg evolves with 
time. This is the most general case and in reality, for k-essence models, it is likely that these two parameters are 
epoch dependent. 

We will study two cases, first with the following evolution of w{a) = —1 + 0.2 (^jj^^, and the two functions 

Ce = (iTa) ^^'^ = (i^)- The choice of w{a) corresponds to the case where it evolves from —1 to —0.9 at 

the present epoch The second case we consider is with the same two functional form of Cg but with w{a) given by 

w{a) — —0.8 — 0.2 ( is initially -0.8 and evolves to —0.9 at the present epoch. The suppression of CDM power 

spectrum as compared to ACDM model in these two cases are shown in the Fig. 7. it is evident from this figure that 
the behavior of suppression with k is similar in the both cases, with suppression being maximum at the smaller scales 
than at larger scales. Also it is clear from the two figures that suppression is more in the case where w is initially 
—0.8 than the case where initial value is -1, reiterating the point that as we move away from value of w = —1 the 
the effect of dark energy perturbations is more important. Finally, we infer that for different functional form of Cg(a) 
for a given evolution of w{a) does not have any considerable effect on the matter power spectrum. In fact it can be 
observed from Fig. 7 that changing the functional form of w{a) for a given evolution of Cg(a) influences the matter 
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FIG. 7: Supression in matter power spectrum as compared to that in ACDM model in the case when both and w are epoch 
dependent. In the left panel we have chosen a form of w{a), which is initially —1 and asymptotically goes to a constant value 
—0.8 and in the right panel w is chosen such that it is initially —0.8 and asymptotically approaches —1. 



power spectrum more severely than the reverse case. 



VI. THE ROLE OF PERTURBATION IN DARK ENERGY 



From the discussions in the preceding section, it follows that for a given evolution w{a), the effect of different 
functional form of on the CDM power spectrum is significant only at the intermediate scales (see Fig.l to Fig. 7). In 
fact, at scales much smaller (or larger) than the Hubble radius, the suppression of CDM power spectrum with respect 
to the ACDM model is nearly independent of the functional form of Cg (a) . The reason for this behavior is discussed 
in this section. 



At the sub-Hubble scales 



At scales much smaller than the Hubble radius, the perturbation in dark energy is negligibly smaller than the 
corresponding perturbation in dark matter [23]. This result follows from the numerical solutions of Eqs.(32) to (34). 
In Fig. 8, we have plotted the ratio (5d/(l + w)dm at the present epoch for different length scales of perturbations. It 
follows from this figure that at scales much smaller than the Hubble radius (fc >> O.OOlh/Mpc), the perturbation in 
dark energy is negligibly smaller in comparison to that in dark matter. In fact, at these scales it can be assumed that 
the dark energy is homogeneously distributed. 

Since we have assumed dark energy to be a scalar field, homogeneous dark energy at these scales reflects to the 
fact that scalar field fiuctuations at these scales are negligibly smaller than other perturbation variables such as Sm, 

etc. With this assumption that the dark energy is homogeneous at scales much smaller than the Hubble radius, it 
turns out that the evolution of the Bardeen potential $ at these scales is governed by the following equation 

$ + 4i7$ + (^2H + 3H^ + 3H^ (^^|^) * = 0- (51) 

The above equation corresponds to the evolution of $ in a system of dark matter and homogeneous scalar field 
dark energy. In the case of ACDM model, the corresponding equation for $ is given by 



$A + 4iJ<i>A + (2H + 3ij2 ) <I>A = 0. 



(52) 
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FIG. 8: The raio of perturbation in dark energy to the perturbation in dark matter scaled by a factor 1/(1 + w) is plotted 
in this figure. It can be observed that at length scales much larger than the Hubble radius this ratio approaches unity. The 
reason for this is that at these scales the incurvature perturbation vanishes. Therefore, at these scales = (1 + 'w)Sm- 



This evolution equation for $ in ACDM model is valid at all length scales. However, it should be noted that Eq.(51) is 
strictly valid only at scales much smaller than the Hubble radius where dark energy can be assumed to be homogeneous. 
At these scales, the perturbation in dark matter is given by 

1 \ e 



^™ = - ( TZTT^ ) -2 * (53) 



Therefore, the percentage suppression in CDM power spectrum due to dark energy in comparison to ACDM model 
(A% defined in Eq.(48)) would become 

$2 _ ^2 

A% = ^ X 100 (54) 

From the numerical solution of Eqs.(51) and (52), it follows that for w = —0.8, the value of A% is —13.39 whereas 
for w — —0.9, it turns out to be —6.17. These values matches with the one shown in Fig.l at k >> O.OOl/i/Mpc, 
which was plotted without assuming a priori the homogeneity of dark energy at these scales. 

The evolution of Bardeen potential $ at scales much smaller than the Hubble radius is independent of Cg (see 
Eq.(51)). The corresponding growth of CDM perturbation at these scales given by Eq.(53) is also independent of 
Cg. Therefore, at these scales the suppression in CDM power spectrum is independent of the dark energy speed of 
sound Cg. As emphasized earlier, the primary reason for this is the fact that at these scales dark energy is nearly 
homogeneous. It is only through the perturbation in dark energy that Cg plays a role in suppressing the CDM power 
spectrum compared to that in ACDM model. 



At the super-Hubble scales 



At scales comparable to Hubble radius and beyond (fc < O.OOl/i/A/pc), the perturbation in dark energy cannot, in 
general, be neglected in comparison to that in dark matter. This evident from Fig. 8, which shows that in the limit 
k « O.OOl/i/Mpc, the perturbation in dark energy becomes 

5d = {l + w)5m (55) 

If w ^ —1, then perturbation in dark energy can become comparable to that in matter at these scales [k < 
O.OQlh/Mpc). Eq.(55) also follows from the fact that at these scales, the isocurvature perturbation or the total 
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non-adiabatic pressure perturbation vanishes. In addition to this, the intrinsic entropy perturbation of dark energy 
also vanishes at these scales. The intrinsic entropy perturbation of dark energy is proportional to its non-adiabatic 
pressure perturbation defined as 

^P„ad = ^Pd - clSpd (56) 

For scalar field dark energy, it <5p„„^ can also be expressed as 

^Pna, = (Ce - cl) [5d - 3H^1 + w)a'ud\ pM- (57) 

At the super Hubble scales, the Sp^^j^ of dark energy vanishes not because of the fact that Cg = but because 

5d = 3H^{1 + w)a^Ud (58) 

In fact, it can be verified that in the limit fc — >■ 0, Eqs.(55) and (58) are consistent with the covariant conservation 
equation for dark matter and dark energy. 

With the fact that the non-adiabatic pressure perturbation of dark energy vanishes at scales much larger than 
the Hubble radius, it turns out that the evolution of Bardeen potential at these scales is governed by the following 
equation 

^+(^ + 3cl{l + w)( -JM^L—) ) + (2H + 3H' + Sif^c^ (1 + w) ( -J^^^^L—] U = 0. (59) 

The evolution of $ in ACDM model follows from the above equation if we substitute w = —1. In the limit k << 
O.OOlh/Mpc, neglecting the k"^ term in time-time linearized Einstein's equation {6G° q = SttGST^q,) together with 
the fact that Sa = {1 + w)Sm, turns out that the equation for dm becomes 

\1 + w{a)ndia) J y H J 

The percentage suppression in CDM power spectrum can now be evaluated by numerically solving Eq.(59) and 
then evaluating Sm using Eq.(60). It turns out A% defined in Eq.(48) for w = —0.8 at these scales {k 0) is —7.7. 

and for w = —0.9 it becomes —3.7. These values are consistent with those in Fig.l for k << O.OOlh/Mpc where we 
have solved the exact perturbation equation without imposing any assumption. 

Hence, at scales much higher than the Hubble radius, the effective speed of sound of dark energy does not influence 
the CDM perturbations. This is because the intrinsic entropy perturbation of dark energy vanishes at these scales. It 
is the intrinsic entropy perturbation of dark energy which carries the information of the ESS parameter Cg(a). Thus, 
it is only at the intermediate scales (scales around k ^ O.Olh/Mpc) where the effect of different evolution of Cg(a) is 
more pronounced on the suppression of CDM power spectrum with respect to ACDM model. 



VII. SUMMARY AND CONCLUSIONS 



In this paper we have investigated the influence of perturbation in dark energy on the matter power spectrum in 
models where the effective speed of sound Cg(a) of dark energy evolves with time. We first presented a method of 
reconstructing the lagrangian density of dark energy of the form £(X, ^) = F{X) — F(<?f)) from a given evolution of 
Cg(a) and the equation of state parameter w{a). This illustrates the fact that these scalar field dark energy models, 
in principle, allows a wide class of solutions with different functional form of w{a) and Cg(a). In fact, in these models 
evolution in Cg(a) can be independent of that in w(a). 

We then investigated the growth of CDM perturbation as influenced by the perturbations in dark energy in models 
where (i) both Cg and w are constants, {ii) either of them are epoch dependent and (Hi) both of them are epoch 
dependent. It is shown that in all these cases, the CDM power spectrum is generically suppressed in comparison to 
that in ACDM model. The degree of suppression at different length scales does in fact depend on the behavior of 
w{a) and Cg(a). 

When the equation of state parameter of dark energy is constant, it is shown that the percentage suppression in 
CDM power spectrum with respect to the ACDM model decreases with increasing length scale of perturbation (see 
Fig.l). However, this is not generically true for any evolution of w(a). In fact, it is found that for slow evolution of 
w{a) from —1 to a value close to it (with w + 1 > 0), the percentage suppression in CDM power spectrum can increase 
with increase in the length scale of perturbation (see, left panel of Fig. 4). 
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Primarily, in this paper wc have compared the percentage suppression in CDM power spectrum in dark energy 
models where its effective speed of sound Cg(a) increases with scale factor with those where it decreases. It is shown 
that the effect of different evolution of Cg(a) of dark energy on the matter power spectrum for a given evolution of 
equation of state parameter w(a) is not as much significant as compared to the reverse case, viz. different evolution 
of w(a) for a given evolution of Cg(a) (as shown in Fig. 5 to Fig. 7). This illustrates the fact that the effect of equation 
of state parameter w{a) of dark energy on the CDM power spectrum is much more severe than its effective speed of 
sound Cg(a). 

Further, it is also shown that the effect of different evolution of Cg(a) for a given evolution of w{a), on the suppression 
of CDM power spectrum is more pronounced only at the intermediate scales at around k ^ O.Olh/Mpc. In fact it is 
observed that the degree of suppression of CDM power spectrum with respect to ACDM model is nearly independent 
of Cg(a) at scales much smaller and larger than the Hubble radius. The reason for this behavior at scales much 
smaller than the Hubble radius is as follows. At these scales, it turns out that the perturbation in dark energy is 
negligibly smaller than the corresponding perturbation in matter. One can effectively approximate dark energy to be 
homogeneous at these scales. Since it is the perturbation in dark energy which carries the information of its effective 
speed of sound, the suppression of CDM power spectrum at these scales will be independent of any functional form 
of Cg(a). The suppression at these scales will primarily be governed by the functional form of w(a) of dark energy. 

Even at scales much larger than the Hubble radius, suppression of CDM power spectrum with respect to ACDM 
model is nearly independent of Cg(a). At these scales the perturbation in dark energy cannot be neglected in compar- 
ison to matter perturbation if w{a) deviates from —1 (as shown in Fig. 8). However, it is found that at these scales 
both the total entropy perturbation (or the isocurvature perturbation) and the intrinsic entropy perturbation of dark 
energy vanishes. It is in fact the intrinsic entropy perturbation of dark energy which encodes the information of its 
effective speed of sound. Therefore, at these scales, although there is dark energy perturbation, it is purely adiabatic 
and consequently its effect on the suppression of CDM power spectrum will be same for different Cg(a). 

In summary, it is illustrated in this paper that the suppression of CDM power spectrum with respect to ACDM 
model both at scales much larger and smaller than the Hubble radius only depends on the form of w{a) of dark 
energy. It is only at the intermediate scales, a non zero value of the effective speed of sound of dark energy leaves 
its imprint on the CDM power spectrum. Precisely determining the effective speed of sound of dark energy, although 
observationally challenging [33], is necessary to understand its nature. We hope that future observations will shed 
more light on the 'darkness' of the dark energy. 
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